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We have conducted an LTE abundance analysis for SAO 40039, a warm post- 
P^ ', AGB star whose spectrum is known to show surprisingly strong He I lines for its 

C/^ ! effective temperature and has been suspected of being H-deficient and He-rich. 

^ ! High-resolution optical spectra are analyzed using a family of model atmospheres 

with different He/H ratios. Atmospheric parameters are estimated from the 
ionization equilibrium set by neutral and singly ionized species of Fe and Mg, 
"^ ! the excitation of Fe I and Fe ll lines, and the wings of the Paschen lines. On the 

assumption that the He I lines are of photospheric and not chromospheric origin. 
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a He/H ratio of approximately unity is found by imposing the condition that the 
^ ! adopted He/H ratio of the model atmosphere must equal the ratio derived from 

the observed He I triplet lines at 5876, 4471 and 4713A, and singlet lines at 4922 
' and 50 15 A. Using the model with the best-fitting atmospheric parameters for this 

j^^ , He/H ratio, SAO 40039 is confirmed to exhibit mild dust-gas depletion, i.e., the 

O I star has an atmosphere deficient in elements of high condensation temperature. 

The star appears to be moderately metal-deficient with [Fe/H]=— 0.4 dex. But 

the star's intrinsic metallicty as estimated from Na, S and Zn, elements of a low 
K/i ; condensation temperature, is [Fe/H]cU=~ —0.2. The star is enriched in N and 

H [ perhaps O too, changes reflecting the star's AGB past and the event that led to 

He enrichment. 
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Introduction 



SAO 40039, also known as IRAS 05040+4820 and BD +48°1220, is judged by several 
criteria to be a post-AGB star: its spectral type (A4Ia), its position in the IRAS color-color 
diagram and the double-peaked spectral energy distribution indicating a detached youthful 
cold dust shell ar ound the central star. Using VBRIJHK photometry and adopting a simple 
dust shell model, iFujii et al.l (120021 ) find a dynamical age of 2460 years for the dust shell with 
dust at 97K, and estimate the star to have a core mass of O.SSM© indicating a low-mass 
progenitor for this PAGB star. 



Klochkova et al.l (120071 ) have presented an abundance analysis of SAO 40039 based on 



high-resolution spectra with a wavelength coverage of 4500-6760A. They have also discussed 
the star's radial velocity changes and also time-dependent differential velocity shifts between 
different classes of absorption lines. Their spectra showed a variable Ha line with two 
emission components. Variable emission was also reported in II/3 and in some lines of Si ll, 
Fe I and Fe ll. But, in particular, these authors pointed out the abnormal strength of the He I 
5876A absorption feature in this star with the 'low' effective temperature of about 8000 K. 
Their abundance analyses conducted with model atmospheres computed for a normal He/H 
ratio (=0.1) gave the He/H ratio of 0.7. 

Detection of He I lines in absorption in seve ral A and F-type PAG B stars has been 
reported previously. The 4471A line was found by IWaelkens et al.l ( 19921 ) in HD 44179, the 
central star of the Red Rectangle. HD 187885, a PAGB star w ith a high C/0 ratio an d 
an enrichment o f s-pro cess elements, shows the He I 5876 A line ( IVan Winckel et al.lll996l ). 
Klochkova et al.l (120021 ) report five He I lines in HD 331319. Detection of He I lines in A and 
F-type supergiants suggest either a He enrichment, if the lines are of photospheric origin, or 
a contribution from a thick hot chromosphere of presumably a normal He abundance. The 
above studies adopt the assumption of a photospheric origin and report the photosphere to be 
moderately H-poor and He-rich. However, these analyses use model atmospheres computed 
for a normal He/H ratio and do not iterate on the construction of the model atmosphere 
and the abundance analysis to obtain agreement between the input and output He/H ratio. 

Our goal in making a fresh analysis of SAO 40039 is to iterate to find the self-consistent 
He/H ratio. In undertaking the analysis, we also adopt the assumption that the He I lines 
originate in a He-rich photosphere and are not chromospheric in origin. 
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Observations and data reduction 



Spectra of SAO 40039 were obtained on ttie niglits of 2007 November 5 and 2011 Febru- 
ary 21 at the W.J. McDo nald Observator y witli tlie 2.7m Harlan J. Smith telescope and the 
TuU coude spectrograph (JTuU et al.lll995l ). The spectra taken at these epochs were found to 
be similar. The spectra correspond to a resolving power of 60,000 (5 kms s~^). Reduction 
of the raw data was performed with the Image Reduction and Analysis Facility (IRAFp 
software package. 

An additional spectrum was obtained on 2011 January 28 using the echelle spectrometer 
at Vainu Bappu Obs ervatory in Kava lur, India giving a resolution of about 28,000 (10 kms 
s~^) in slitless mode ( iRao et al.ll2005l ). Since the projected rotational velocity of the 
star is about 15 kms s^, the lower resolution of the VBO spectrum does not 
affect the stellar absorption profiles and hence our measured equivalent widths 
and the calculated abundances. 



3. Abundance Analysis 



Our spectra resemble closely those described by iKlochkova et al.l ( 120071 ) with regards 
to emission in Balmer lines, the complex structure of the Na D lines and variable profiles 
of the metallic lines, and the equivalent widths of lines unaffected by obvious emission. For 
our analysis we have used only clean, unblended and symmetric absorption lines. 

We have compared th e equivalent widths of lines common between our spectra and those 
of IKlochkova et al.l ( l2007l ). A small systematic difference possibly caused by differences in 
resolution is discernible, but it appears that there were no major variations in the atmospheric 
parameters at these epochs. This is also supported by the fact that the avera ge radial velocity 
of line s in our spectra (-12 kms s~^) lies within the range as reported by IKlochkova et al. 
( 120071 ) in their spectra (-7 to -15 kms s~^). 



Abundance analysis was done first on the assumption that the atmosphere has a normal 
He/H ratio. Model atmospheres for normal He/H ratio (He/H=0.1) were taken from the Ku- 
rucz database!. Atm ospheres of different He/H ratios were computed by the code STERNE 
( iJeffery et al.ll200ll ). The LTE spectrum synthesis code MOOG (2009 version) by ISneden 



^The IRAF software is distributed by the National Optical Astronomy Observatories under contract with 
the National Science Foundation 

•^ http://kurucz.harvard.edu/grids. html 
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( 119741 ) was used with the grid of model atmo spheres taken from the Kurucz database as men- 



tioned above. The LTE code SPECTRUM flJefferv et aLll200l[ ) was used with the STERNE 
models 111 

The procedure begins with the determination of the atmospheric parameters from the 
spectrum, continues with the abundance analysis for He and other elements and ends when 
the model of a certain He/H ratio reproduces the observed He I lines. 



3.1. Method for determining the stellar parameters 

First, we estimated the microturbulent velocity {^t) using Fell lines with small range 
in lower excitation potentials (LEP) (2.6—3.2 eV) as these lines have a good range in their 
equivalent widths thereby reducing any temperature dependence of the estimated ^f The ^t 
is found from the standard requirement that the abundance be independent of the measured 
equivalent width. The effective temperature (Tcs) was estimated by the requirement that 
the abundance of a given species be independent of a hue's LEP. This step was conducted 
independently for Fel and Fell lines as they are both well represented in the spectrum 
and show a range in their LEP's. Finally, the surface gravity {logg) was estimated by the 
condition that there be ionization balance between the neutral and the singly ionized Fe 
lines. This condition defines a locus in the (Teg, logg) plane and the T^s derived from Fel 
and Fe ll breaks the degeneracy. This exercise is repeated for grids of model atmospheres 
with He:H ratios of 10:90, 30:70 and 50:50. 

A microturbulence ^t = 4.8 ± 1.0 km s^^ is found for all models with acceptable effec- 
tive temperature and surface gravity and the value is insensitive to the He/H ratio of the 
models. The effective temperature T^s is found to be 8000±300K from 26 Fe I and 20 Fe ll 
lines. Within the 300K uncertainty, the temperature is independent of surface gravity over 
a considerable range and is not sensitive to the He/H ratio of the model atmosphere. 

The Teff, logg locus found from the ionization balance for Fe is illustrated in Figure 
1 for the model atmosphere grid with He:H = 50:50. With Teg = 8000±300 K from the 
excitation of Fel and Fell, the log (?=0.75±0.25 cgs is found. The surface gravity changes 
slightly with He:H running from 0.94 for He:H=10:90 through 0.83 for He:H=30:70 to the 
above result for He:H=50:50. A second locus is provided by Mgl and Mgll lines - see upper 



"^Note that the steUar parameters and the abundances derived by adopting the Kurucz models and the 
code MOOG for He/H = 0.1 are in excellent agreement with results from the STERNE models and the LTE 
code SPECTRUM, also for He/H = 0.1. 
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panel of Figure 1. The Mg and Fe loci are in good agreement. A third locus is offered by 
the Paschen lines (see below). 



3.2. The He/H ratio 

Our spectra confirm th e presence of the 03 587 6A He I line with an equivalent width very 



similar to that reported by iKlochkova et al.l (120071 ). The strength of this line suggests that 



several other lines should be present unless the excitation of the D3 line is highly peculiar. 
Our search also provided the additional triplets (4471 and 4713A) and two singlet lines (4921 
and 5015A). Since the He I at 4713A is blended with the Fell line, its contribution has been 
included while synthesizing the He I profile at 4713A as shown in the Figure 2, we get a 
consistent He/H ratio of approximately unity. Other possible He I lines are predicted to be 
below the detection limit, or blended with other lines or inaccessible due to inter-order gaps. 

The He I lines were computed for model atmospheres with {Tcs, logs', 6) = (8000, 0.75, 
4.8) and with the appropriate model for the He/H ratios: 0.1 (normal; solar abundance), 
0.4 (30/70), and 1.0 (50/50). The synthesized profiles are convolved with the instrumental 
and the stellar rotation profiles, before matching with the observed spec trum as described i n 



Pandeyl (120061 ). The data for computing the He I profiles were taken from lJeffery et al.l (J200ll ). 
A projected rotational velocity of about 13—16 km s^^ is estimated by using unblended 
moderately strong lines. As shown in Figure 2, the best fits to the observed He I profiles 
at 5875, 4471, 4713, 4922 and 5015A were obtained for He/H ratio of 1.0 (50/50) with an 
uncertainty of approximately ±0.2. 

To verify the model atmosphere parameters, we explore fits to the Paschen lines. The 
Balmer lines are affected by emission and blends with other lines for their line profiles to 
be useful as monitors of the atmospheric parameters but the Paschen lines appear free of 
emis sion and blends. The Stark broadening tables for the Paschen lines of hydrogen are taken 



from iLemkd (119971 ). In the lower panel of the Figure 1, we show observed and synthesized 
Paschen profiles near 8600A for the model with He/H=1.0 and T^g = 8000K and for surface 
gravities logg = 0.5, 0.75 and 1.0. The best fitting profile is for \ogg of 0.75. We also 
synthesized the Paschen lines in the 8400A region for the above mentioned temperature 
and gravities and we get consistent results. These profiles are also sensitive to the adopted 
effective temperature with a lOOOK increase approximately mimicking a 0.5 dex increase in 
surface gravity. The Paschen lines give the locus shown in the upper panel of Figure 1. As the 
Paschen limit is approached, the Paschen lines overlap and depress the local continuum. In 
addition, the location of the continuum is rendered more challenging because of the echelle's 
blaze function. Yet, the synthesized profiles closer to the limit are able to reproduce the 
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observed profiles quite well - see lower panel of Figure 1. 

Our demonstration that SAO 40039 is H-deficient and He-rich rests on the assumption 
that the He I absorption lines originate in the stellar photosphere and are not a product of a 
stellar chromosphere. (A second assumption is that the real stellar photosphere approximates 
the model atmospheres computed by classical procedures). A chromospheric origin is not 
trivially dismissed. One rec alls that the D3 line appears in absorption in sp ectra of F dwarfs 



with active chromospheres ( IWolff fc Heasleyl Il984j : iDanks fc LambertI 119851 ) ; the correlation 



between D3 equivalent width and X-ray luminosity points strongly to a chromospheric origin. 
A chromospheric origin for He absorption is required to account for the presence of He 1 10830 
A in spec tra of warm-cool giants ari d supergiants where it may appear in absorption and/or 



emission (jO'Brien 8z LambertI Il986l ) and cannot be of photospheric origin. 



The possibility of the photospheric origin of the He I lines in our spectrum is supported 
by several facts. Firstly, the lines have the photospheric velocity. Secondly and more im- 
portantly, the equivalent widths of all He I lines - singlets and triplets - are reproduced 
by a photospheric model. If the lines originated in the chromosphere, one would expect 
a difference between singlets and triplets, as previously noted by Klochkova et al. (2002) 
in their discussion of the He I absorption lines found in HD 331319, a star similar to SAO 
40039. While our analysis do not preclude a chromospheric origin, we agree with Klochkova 
et al. who remark 'the hypothesis of a photospheric origin for the He I lines in HD 331319 is 
not rejected'. Certainly, there must be a suspicion that SAO 40039's atmosphere is mildly 
H-deficient such that He/H= 1 now, a value indicating appreciable loss of H and enrichment 
of He since the star was a main sequence star with He/H ~ 0.1. 



4. Discussion and Results 

Using the adopted model parameters (Teg, logg, ,^t) = (8000, 0.75, 4.8) with He/H of 1.0, 
the measured equivalent widths of different species were used to derive the abundances. In 
Table 1 we present the abundances derived with the SPECTRUM code using the He/H=1.0 
model and the abundances generated with the code MOOG using the He/H=0.1 model. 

A detailed line list (see for SAMPLE Table 2, that lists some hues of SAO 40039) used 
in our analysis has been presented electronically. The line list provides the lower excitation 
potential (x) for each line , the log gf value, sources of loggf, the measured equivalent 
widths {W\) in mA the abundance (log e) derived from each line for the adopted model 
atmosphere. 



The solar abundances for all the elements have been taken from lAsplund et al.l (120091 ) 
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SAO 40039 seems to be mildly iron-poor. Klochkova et al.'s analysis assuming He/H= 
0.1 gave a solar Fe abundance, a value essentially confirmed by our analysis for the model 
atmosphere also with He/H = 0.1. A principal effect of the higher He abundance is to 
reduce the continuous opacity per gram and so demand a lower iron abundance to match 
the same observed line strength. Inspection of the [X/Fe] entries in Table 1 points to several 
anomalies when the composition of an unevolved disk star is taken as the reference. The 
two most extreme anomalies are nitrogen which is very overabundant ([N/Fe] = +1.1), and 
aluminium which is markedly under- abundant ([Al/Fe] = —0.7). Closer study shows that 
sodium may be overabundant ([Na/Fe] = +0.3) and a suite of elements are mildly under- 
abundant: Ca, Sc, Ti, Zr and Ba with [El/Fe] ^ —0.35. 

This abundance pattern is largely reminiscent of the pattern exhibited by stars af- 
fected by dust-gas separation, i.e., the photosphere is dominated by accretion of gas but 
not dust from a cool envelope, possibly a circumbinary dusty disk. Dust-gas separatio n 
is most remarkably exhibited by post-AGB binaries such as HR 4049 ( IVan Winckell l2003l ) . 
In such cases, the abundance anomalies are correlated with the condensation temperature 
(Tc): elements with high T^ are most under- abundant. In Figure 3, we plot [X] (where [X] 
=loge(X)*— loge(X)0) versus the Tc computed by Lodders (2003) for a solar composition 
mixture. Elements N and O depart from a general tendency for [X] to decrease with in- 
creasing Tc. In presenting Figure 3, the adopted reference has been that [X/Fe] ~ 0.0 for all 
elements in an unevolved star with [Fe/H] ~ —0.4. This is essentially true f or a thin disk star, 



althou gh a dispersion in heavy elements at a given [Fe/H] was noted by lEdvardsson et al. 
( 119931 ). If SAO 40039 is a thick disk star, the a-elements in Figure 3 should be plotted 
with [X] about —0.2 dex less than shown in the figure but the overall trend and the scatter 
about that trend will be unaffected. The Na, S, and Zn abundances suggest that the initial 
metallicity of the star was [Fe/H]o — —0.2. 

The high N abundance is, probably, a consequence of extensive burning of H to He 
by the CN-cycle with concomitant conversion of C to N. Survival of C at about its initial 
abundance suggests that N synthesis occurred primarily after the atmosphere was enriched 
with C at the third dredge-u p. One possibility is H-consumption during a final He-shell flash 
( lHerwig||2001t lBlockerll200l[ ). The O is overabundant with respect to the initial abundance, 
possibly it was enriched by the third dredge-up and temperatures during the late episode of 
H-burning were too low to destroy O or non-LTE line formati on can reduc e the C, N and 
O abundances. For an effective temperature of about 8000K, IVennI ( 1l995l ) give non-LTE 
correction to C and N abundances from Cl and Nl lines of about —0.4 dex for normal 



supergiants. iTakeda fc Takada-Hidail (119981 ) have estimated NLTE corrections for O I to 
be —0.1 dex at this effective temperature for supergiants. After applying these non-LTE 
corrections, the above sketch of a possible origin for SAO 40039 is little affected. 



5. Concluding remarks 



Though, the existence of extremely H-deficient st ars with He/H ratio of about 10,000, 
showing He I hues in their absorption spectra, is known ( iPandey et al.ll200ll : iPandey &: Reddy 
20061 ) ■ SAO 40039, from our estimated He/H ratio of about unity, appears to be the first 
detection of a mildly H-deficient star. iLambertl (119961 ) had commented on the possibilities 
of such stars and on the difficulty in detection for stars too cool to show He I lines. With an 
effective temperature of about 8000K, SAO 40039 seems to be responding to this call. 

It is very likely that other post-AGB stars may belong to the mildly H-deficient class. 
Several examples like SAO 40039 have detectable He I lines. The question then arises; are 
post-AGB stars of a temperature too low to provide He I lines also He-rich and how can 
their H-deficiency be unmasked? Exploration of post-AGB stars showing He I lines need to 
be continued with model s of appropriate H/He rat ios. Additionally, non-LTE abundance 
analyses as described by iPandey &: LambertI ( 120111 ) should be performed to further refine 
knowledge of the chemical compositions. 

We are thankful to the VBO observing support team, in particular to Mr G. Selvakumar, 
for obtaining a spectrum of SAO 40039 using the VBT echelle spectrometer. DLL thanks 
the Robert A. Welch Foundation of Houston, Texas for support through grant F-634. 
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Table 1. Comparison between abundances estimated with MOOG using He/H=0.1 model 
and the abundances computed with SPECTRUM using the hydrogen deficient (He/H=1.0) 

model 
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6.21 


-0.22±0.00 


1 


5.64 


-0.27±0.00 


1 


0.4 


-0.36±0.00 


1 


Sell 


3.15 


2.89 


-0.35±0.17 


5 


2.42 


-0.30±0.15 


5 


0.3 


-0.05±0.16 


2 


Till 


4.95 


4.83 


-0.21±0.23 


20 


4.32 


-0.20±0.20 


20 


0.3 


-0.24±0.16 


12 


VII 


3.93 


4.21 


+0.19±0.04 


2 


3.74 


+0.24±0.03 


2 


0.2 






CrI 


5.64 


5.62 


-0.11±0.13 


3 


5.11 


-O.lOiO.lO 


3 


0.4 


-0.06±0.23 


2 


Cr II 


5.64 


5.83 


+0.10±0.11 


9 


5.33 


+0.12±0.10 


9 


0.2 


+0.09±0.06 


4 


Mnll 


5.43 


5.73 


+0.21±0.00 


1 


5.24 


+0.24±0.00 


1 








Fel 


7.50 


7.54 


-0.05±0.18 


26 


7.02 


-0.05±0.17 


26 


0.4 


-0.01±0.19 


7 


Fell 


7.50 


7.65 


+0.06±0.19 


20 


7.13 


+0.06±0.17 


20 


0.2 


+0.01±0.14 


9 


Nil 


6.22 


6.16 


-0.15±0.15 


4 


5.63 


-0.16±0.14 


4 


0.3 






Znl 


4.56 


4.81 


+0.16±0.30 


2 


4.28 


+0.15±0.23 


2 


0.4 






SrU 


2.87 


3.01 


+0.05±0.07 


2 


2.38 


-0.06±0.03 


2 








Zr II 


2.58 


2.13 


-0.54±0.08 


2 


1.65 


-0.50±0.04 


2 


0.3 






Ban 


2.18 


1.85 


-0.42±0.15 


2 


1.30 


-0.45±0.11 


2 




-0.35±0.00 


1 


EuII 


0.52 


0.40 


<-0.10±0.00 


1 


0.15 


<+0.06±0.00 


1 









^MOOG abundances computed using the He/H=0.1 model for the McDonald spectra. 

''SPECTRUM abundances computed using the He/H=1.0 model for the McDonald spectra. 

'^SPECTRUM abundances computed using the He/H=1.0 model for the VBO spectra. 

'^A corresponds to the square root of the sum of the squares of the abundance errors due to uncertainties in the stellar 
parameters: A Teg, Alogg and A^t 
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Table 2. Detailed line list for SAO 40039 ^ 



Line 


X 
(eV) 


logc?/ 


Wx 
(mA) 


loge 


He I A5875.61 


20.96 


+0.739^^ 


Synth 


11.45 


He I A4471.47 


20.96 


+0.053*^ 


Synth 


11.40 


HelA4921.93 


21.21 


-0.435^ 


Synth 


11.45 


HelA5015.67 


20.61 


-0.819^ 


Synth 


11.50 


HelA4713.13 


20.96 


-0.975^ 


Synth 


11.40 


Ci A4775.87 


7.49 


-2.304^ 


35 


8.41 


Cl A5039.05 


7.48 


-2.100^ 


33 


8.15 


Cl A4770.02 


7.48 


-2.437^ 


30 


8.45 


Ci A4766.67 


7.48 


-2.400^ 


26 


8.34 


Cl A4826.80 


7.49 


-2.140^ 


24 


8.04 



'^NIST database (http://physics.nist.gov/PhysRefData/ASD/lines_forni.htnil ) 

^Table 2 is published in its entirety in the electronic edi- 
tion of ApJ. A portion is shown here for guidance regarding 
its form and content. 



